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Abstract 

We solve the Faddeev equations for the two meson-one baryon system TrvriV and coupled 
channels using the experimental two-body t-matrices for the ttN interaction as input and 
unitary chiral dynamics to describe the interaction between the rest of coupled channels. In 
addition to the A^*(1710) obtained before with the vrvrA^ channel, we obtain, for J'" = 1/2+ 
and total isospin of the three-body system / = 1/2, a resonance peak whose mass is around 
2080 MeV and width of 54 MeV, while for / = 3/2 we find a peak around 2126 MeV and 42 
MeV of width. These two resonances can be identified with the A^*(2100) and the A(1910), 
respectively. We obtain another peak in the isospin 1/2 configuration, around 1920 MeV 
which can be interpreted as a resonance in the A^ao(980) and V/o(980) systems. 

1 Introduction 



Recent developments around three-body systems with two mesons and one baryon using chiral 
dynamics have brought new light into the nature of the = 1/2"'" baryonic resonances. The 
study of such systems with strangeness 5* = — 1 produced resonant states which could be identified 
with the existing low lying baryonic = 1/2'^ resonances, two A and four S states [H E], [3]. 

* aniartine@ific.uv.es, kanchan@teor.fis.uc.pt, oset@ific.uv.es 
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Similarly, in the case of the S = sector the A^*(1710) appears neatly as a resonance of the tttcN 
system, with or without including its coupled channels within SU(3) [1]. Developments along the 
same direction produced a resonant state of (f)KK [5] which could be identified with the X(2175) 
resonance reported at BABAR [6l [7] and later on at BES [8] . The study of the three-body systems 
was done using Faddeev equations (FE) in the coupled channel approach. While most conventional 
studies of three-body systems use potentials in coordinate space, usually separable potentials to 
make the solution of the FE feasible, the approach of [H [2] used two particle amplitudes generated 
within the unitary chiral approach in momentum space. Yet, the most novel finding in these 
works was the realization that, for s-waves and in the SU(3) limit, there was an exact cancellation 
between the off shell part of the two-body amplitudes and the three-body forces generated by 
the same chiral Lagrangians. To be more precise, the on shell amplitude means that the s-wave 
amplitude is calculated as a function of the Mandelstam variable s imposing = m? for the 
external momenta of the two body amplitudes. When these lines are inside the Faddeev diagrams 
where some line can be off shell, the full amplitude is separated into this "on shell" part plus 
and "off shell" part which goes as — m? for mesons and g° — E{q) for baryons and vanishes 
when the external lines are on shell. This off shell part contains an inverse particle propagator 
and cancels one particle propagator rendering a Faddeev diagram with two two-body t-matrices 
into a three-body contact term, which has the same topology as genuine three-body interactions 
that stem from the chiral Lagrangians and cancel them exactly. As a consequence, one needs 
only the on shell two-body t-matrices and can ignore these three-body forces. This finding is 
novel for such studies and simplifies the work technically, although not much, since loops involve 
a changing s-variable, and consequently the s-dependent t-matrices must be inserted into the loop 
functions. This makes this approach different and technically more involved than the study of 
the two body interaction, where using arguments of the N/D method one can factorize on shell 
amplitudes outside the loop functions which involve only two hadron propagators [9], [10]. The 
strongest value of that finding in the three-body problem is that the results do not depend upon 
the off shell extrapolations of the amplitudes which is a source of uncertainty in the three-body 
calculations that rely upon a potential. Indeed, it is well known that given a certain physical 
amplitude, on shell by nature, one has an infinite number of potentials that give this amplitude 
upon solving the Schrodinger equation. The differences between the different potentials will only 
show in the off shell extrapolation of the amplitudes. However, this information enters the solution 
of the Faddeev equations and, hence, different potentials leading to the same on shell amplitude 
will provide different results upon solution of the Faddeev equations. 

The problem stated above is most probably the main reason why recent works dealing with the 
KNN system lead to quite different results in the binding and the width. In this sense, we find 
a series of works based on Faddeev equations which lead to relatively large binding, of the order 
of 50 — 70 MeV [TTl [T2| [T3l [H], while other works based on variational methods lead to smaller 
bindings of the order of 20-30 MeV [13 [ISl [17] . The widths also vary from 50 - 100 MeV. 
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The arbitrariness of the off shell amplitude is also well known in field theory, where the imple- 
mentation of unitary transformations of the fields in the Lagrangian maintains the same on shell 
amplitudes but changes their off shell extrapolation. In this sense it is interesting to note that, 
although the off shell versus three-body cancellation discussed here is not explicitly shown in other 
three-body works using also chiral dynamics |T8l[T9], the approach is invariant upon these transfor- 
mations, indicating that the mentioned cancellations apparently occur in the full calculation [20] . 
A similar independence on the off shell extrapolation has been shown in different reactions like the 
ttN ttttN reaction [21] and the study of the interacting two pion exchange in the NN interaction 
[22]. However, the explicit realization of the off shell versus three-body forces indicates that one 
can neglect the three-body forces from the beginning, certainly simplifying the approach, and use 
only the two-body on shell amplitudes. Even more, these on shell amplitudes can be obtained from 
experiment and one can omit having to do a theory for the two-body interaction. There is a small 
caveat there, since sometimes in the loops one will need "on shell" amplitudes below threshold. 
This looks like a contradiction, but we made it clear the meaning on the on shell amplitude needed 
in the Faddeev equations, which is the one where = for the external momenta. Provided one 
has a suitable parameterization of the amplitude, the extrapolation below threshold fulfilling this 
condition is not a difficult task to accomplish. In many cases, like in the present one that we shall 
discuss here, one needs the information well within the physical range and the extrapolation is not 
even needed. This picture presented here is rather novel and the purpose of the present paper is to 
show how it works and how it can help whenever the theoretical models are not accurate enough. 

With the perspective given above we shall tackle here the investigation of three-body systems 
with two mesons and a baryon with strangeness S = 0. The problem was already discussed in 
[1], where the iV*(1710) was found as a resonant state of ttttN. It was also found there that the 
implementation of other coupled channels barely changed the results obtained with the base of 
the nnN states alone. Yet, there are other = 1/2+ states, like the iV*(2100) and the A(1910), 
which do not appear with that base and the use of the amplitudes obtained with the lowest order 
chiral lagrangians. From the work of [23] we know that the chiral unitary approach using the 
lowest order chiral Lagrangian provides a fair amplitude up to ^/s = 1600 MeV but fails beyond 
this energy. For instance, the A^*(1650) does not appear in the approach. As a consequence, any 
three-body states which would choose to cluster a vrA^ subsystem into this resonance would not 
be obtained in the approach of [4j. In the present work we shall give the step to use experimental 
TriV amplitudes and will show that in this case we reproduce the A^*(1710) resonance without 
practically any modification with respect to |@], but the use of a more realistic nN interaction 
at higher energies leads also to the generation of the A^*(2100) and the A(1910) resonances as 
three-body systems of two mesons and one baryon in coupled channels. 
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2 Formalism and Results 



We follow the method developed in [H [21 [Sj HI [5] to calculate the three-body T-matrix and search 
for resonances. In [H [21 [3l HI [5] a coupled channel Bethe-Salpeter equation is solved to calculate 
the required two-body t-matrices with the potentials obtained from chiral Lagrangians. These 
t-matrices, which contain the information of the two-body resonances, are then used as an input 
to solve the Faddeev equations in a coupled channel approach. The Faddeev equations, 

r = f6^{k^ - ki) + fg'^T^ + fg'^T'', j y^k = 1,2,3 (1) 
in our formalism, have been reformulated to 

T = f5^ik;-ki) + ri + (2) 

where T]j satisfy the equations 

tH = fg'H^ + f [G'^^T^j^ + G'^'Tii] . (3) 

Eqs. Q are matrix equations since we solve them for coupled channels. In Eqs. ([3]) , t"'s are 
the two-body t-matrices for the interaction of the pair of particles {ml) , with n ^ m ^ I {n is 
the index of the spectator particle) and the elements of the g'^^ matrix can be written in a general 
form as 

g^^{ki k,) =(f[^] Kr-^ 1^1^3 = 1,2,3 (4) 

V.=i ^Er) ^s- E,{k, ) - Ei{k, + k,) - E^ikj) 



where D is the number of particles propagating between two t-matrices. Following the normal- 
ization of [21], Nr = 1 for a meson and A^,. = 2Mr for a baryon with Mr being the mass of the 

baryon and ki (kj) represent the momentum of the zth (jth) particle in the final (initial) state (in 
the global center of mass system). 

The G^^^ matrix in Eqs. ([3]) is defined in terms of the product of two matrices as 

r d^k" -> 

G^'^ = J j^~9''isim, k")r^'{k", k'^, kk, sfy) i^j,j^k = l,2,3 (5) 

with I ^ m ^ ^ k' ^ j , where the elements of the g^^ matrix are given by 

= ^ , (6) 

2Ei{k") 2Era{k") - Ei{k") - Era{k") + le 
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and the F"^^^ matrix is defined as 

F^^\k", k^, h, st:i) = f{st:i,) X g^\t', h) X y\k'^, h)]-' X [t\s,k')]~'. (7) 

In Eq. ([6]), Ei{k") = \J {k"Y + and ^Jsim is the invariant mass of the [Im) pair (with 
i ^ I ^ m). We solve Eqs. ([3]) as a function of the total energy of the system, -y/i, and the 
invariant mass of the particles 2 and 3, a/s23. The rest of the kinematical variables, i.e., the other 
Sjj's, energies and momenta of the three particles in different reference frames, are defined in terms 
of these two variables as shown in detail in [U |Tj. In Eq. ([7]) the k" superindex on sf/^i., indicates 
that the definition of this invariant mass depends on the loop variable. 

Eq. ([7]) is derived in such a way that diagrams up to three t matrices are calculated exactly by 
multiplying the contribution of diagrams with two t matrices by the loop function of Eq. ([5]) and 
one more t matrix. This same iteration procedure is used to go from diagrams containing three 
to four t matrices and so on. The procedure was proved to be numerically accurate in [1] . As one 
can see, the formalism has been developed in such a way that all the loop variable dependence is 
contained in the (^''""-functions and thus Eqs. ([3]) are an algebraic set of equations. 

We solve Eq. ([3]) for the ttuN system and coupled channels. The present work benefits from 
the previous study of the ttttN system and coupled channels where the dynamical generation of 
the A^*(1710) was found. This state was found when the total isospin 1/2 three-body T-matrix was 
evaluated by adding a nucleon to the two pions interacting in isospin zero in the energy range of 
the cr-resonance. The A^*(1710) was thus interpreted as a resonance in the vrvrA^ system, where the 
two pions rearrange to form the cr-resonance. The total energy range studied in [1] corresponded 
to a variation of the invariant masses of the vrA^ pairs up to ~ 1550 MeV. The calculations in [4| 
were limited to this energy range because the input vrA^ t-matrix used in that work was taken from 
[23] which reproduces the vrA^ scattering data well up to about 1600 MeV. 

The motivation of this work is to extend the calculations made in [1] to higher energies by 
including the A^*(1535) and iV*(1650) in the input vrA^ t-matrix and look for the other three-body 
isospin 1/2 and 3/2 states with J-'^ = 1/2"*" in the ttttA^ system and coupled channels. In order to 
do this, we use the experimental L = phase shifts (5) and inelasticities [rj) [25j for the ttN system 
in isospin 1/2 and 3/2 configurations ( Fig. [H [2]) and calculate from them the vriV amplitudes in 
the isospin base (Fig. [3] ) using the relation 

t' = -^/', 1=1/2,3/2 (8) 

with 

I 2iS' _ 1 

where rj^ is the inelasticity, 6' the phase shift, M is the nucleon mass, E is the vrA^ center of mass 
energy and q is the corresponding momentum. 
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Figure 1: Experimental phase shifts and inelasticity for the ttA^ interaction in isospin 1/2. 
5^'\deg.) 





1200 



1400 1600 
E (MeV) 



1800 



1200 



1400 1600 
E (MeV) 



1800 



Figure 2: Experimental phase shifts and inelasticity for the ttA^ interaction in isospin 3/2. 



We require the input two-body t-matrices in the charge base to solve the Faddeev equations in 
our model. For this we use the relations 



23/2 1^1/2 

3 3 ' 



3 3 
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Figure 3: Experimental t-matrices for the vrA^ interaction in isospin 1/2 and 3/2. 



tn^p—^TT+n 



1^3/2 ^ 2 1/2 

3 3 ' 



^TT n — s-TT n 



p — *7r p7 



tn^p—t-K^p 



(10) 



Using these irN t-matrices as input for Eqs. ([3]), we can extend the model for the tttcN 
interaction of |1] to higher energies where the invariant masses of the irN subsystems can be varied 
around 1650 MeV. 

At this point we would like to discuss the cancellation between the off-shell part of the t-matrices 
in the Faddeev equations and the three-body forces, which justifies the use of the experimental 
amplitudes in our approach. These cancellations have been illustrated in all detail in the appendix 
of |1] for the systems of two mesons and one baryon and in the appendix of [S] for the one vector-two 
pseudoscalar systems. The proof proceeded taking the potentials (tree level amplitudes) derived 
from the lowest order chiral Lagrangians. The extention of the proof made in the appendices of 
[H [5] to the corresponding one using t-matrices is straight forward, since the t-matrices would be 
generated by further iterations of the tree level amplitudes in the Faddeev diagrams, as done in 
P, |23j, where the off-shell part of the potential in these iterations is reabsorbed in constants of 
the on-shell potential. Hence the cancellations are guaranteed when iterations are done to obtain 
Faddeev diagrams in terms of t-matrices rather than potentials. This was discussed in [U [5]. 

One could wonder if such a cancellation would also occur in those cases where the higher order 
terms of the Lagrangian would be necessary. Technically our assertion, that one can use only the 
on-shell amplitudes, is rigorous as long as the amplitudes obtained with the lowest order chiral 
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Lagrangian, upon unitarization, can reproduce the experimental data. This seems to be the case, 
for example, in S=-l systems, for the energy range considered here. Indeed, calculations done 
in [26] using higher order terms of the Lagrangians show that the results obtained by using the 
lowest order Lagrangian fall well within the accepted uncertainties in the model. The situation is 
different for 5 = 0, since, as mentioned in the introduction, the results of the calculations done 
with the lowest order chiral Lagrangian already fail beyond the total energy of 1600 MeV of the 
TiN system. Thus we could formally make no claims in this region about cancellations between the 
off-shell part of the T^-matrices and the three-body forces. However, we also insist on the fact that 
the results cannot depend on the off-shell part of the amplitudes, because these are unphysical. 

The finding of the exact cancellation of the off-shell part with the three-body forces is very useful 
because it implies that the Faddeev equations can be solved using only the physical information, 
that is the on-shell amplitudes. This feature certainly must sustain even when one goes beyond 
that realm where the lowest order chiral Lagrangian reproduces the experimental data. It would 
be interesting to study cancellations similar to those found in [U |5] for the present case by using 
higher order Lagrangians but this is beyond the scope of the present work. 

It should be also said, when using higher order terms in the Lagrangians, that although the 
ehmination of the off-shell (unphysical) part is guaranteed, because the results cannot depend on 
unphysical amplitudes, it is not clear that the cancellation mentioned above would not leave some 
finite remanent part. It is also not guaranteed that, apart from three-body forces originating from 
the chiral Lagrangians, there are no other genuine three-body forces which would remain after 
necessary cancellations of off-shell terms. 

However, let us make the following observation. We used a theory suited to the study of the 
TiT^N system up to ^/s ^ 1750 — 1850 MeV and concluded that one can study the system using 
only on-shell amplitudes, which one can get from experiments. Although proved within a certain 
theory, the conclusion that "one can solve Faddeev equations with experimental amplitudes" is 
not linked to any model. With this in mind we make an ansatz that this conclusion should not 
be linked to the theory used to prove it and should be a characteristic of the dynamics of these 
systems for a wider range of energies than the one where we could establish a proof based on a 
particular theoretical framework. Afterall we will only extend our region of energies up to ^/s ~ 
2200 MeV which is not too far from the energies at which the calculations were made earlier. 
Although certainly it is an ansatz at these higher energies, a posteriori, our assumption that one 
can rely solely upon the on-shell amplitudes in the Faddeev approach gets a strong support from 
the results that we obtain in the present work. 

2.1 Exploring the ttttN system 

We first study the ttttN system with total charge zero considering vr'^vr^n, 7t^tt~p, T[^Ti~n, tt^tt+tt, 
and TT~7c^p as coupled channels. We label them as particle 1, 2 and 3 in the order in which they 
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are written above. We calculate the three-body T}1 matrices (Eqs. ([3])) by using, for the nN 
interaction: (a) experimental amplitudes, i.e., Eq. with the phase shifts and inelasticities 
shown in Figs. [H when the invariant mass of vrA^ system is above its threshold (b) and the 
t-matrix obtained from chiral Lagrangian |23] for those ttN total energies which fall below the 
threshold. For the vrvr interaction we use the t-matrix obtained and studied thoroughly in [27] . 
where the dynamical generation of the (j(600), /o(980) and ao(980) resonances was found and 
the theoretical results for physical observables coincided well with the experimental ones. We 
take proper symmetrized amplitudes into account wherever necessary, for instance, for the vr^vr^ 
subsystem in the ti^ti^N channel. 

As in we obtain the total T^-matrix defined as =J2ij(T}l —fg^H^) (see [1\ for this 
definition). In order to identify the nature of the resulting states, we project the T^-matrix on the 
isospin base. An appropriate base is the one where the states are classified by the total isospin 
of the three particles "/" and the total isospin of a subsystem We thus label the states 

in the isospin base as | /, Igub) ■ Obviously transitions between states with same total isospin but 
different isospin of a subsystem are possible. The peaks in the amplitudes are nevertheless seen 
more clearly for some particular isospin of the subsystem, indicating that the dominant structure 
of the state found in the three-body system has a certain value of the total isospin and that of 
the isospin of a subsystem. We can thus write our T^-matrix in the isospin base, in general, as 

{I, huh I T^iV^, y/S23) I I, hub)- 

In order to be consistent with our previous work [4j, we first check if we find an evidence for 
the A^*(1710). For this, we obtain the amplitude for total isospin of the three particles 7=1/2 
and for the isospin of the vrvr subsystem, denoted by J^tt, being equal to zero, i.e., we calculate 
(/ = 1/2, Jtt^ = 0|TR(A/i, ^/sq^)\I = 1/2, Jtt^ = 0) (see Eqs. (30) of [1] for a more detailed definition 
of this isospin base) . 

We find exactly the same peak at 1704 MeV in the squared amplitude as obtained in |4j. In 
this way, we ensure that we reproduce our previous results by using the experimental data for the 
ttN interaction above the vrA^ threshold. With this assurance, we now look for resonances in the 
higher energy region in same or other isospin configurations. 

We now study the irnN amplitudes for the case in which the isospin of the subsystem of particle 
2 and 3, i.e., pion nucleon (and its coupled channels) is 1/2 in the initial as well as the final state. 
To obtain this amplitude we write the vrvriV states in the isospin base as 



1,0)® I 1,0)® I 1/2,-1/2) 
l,0)®<;i/^ |/.iv = 3/2,J^^ = -l/2) 



- I Ln - 1/2, I^N 



/ = 5/2,/^;v = 3/2) 



^/f|/ = 3/2,/.. = 3/2)-^|/ 



(11^ 

-1/2)^ 
1/2,1^^ = 3/2) + 
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Figure 4: The TV* (2100) in the ttttN system with five coupled channels. 



+ ^ I 7 = 3/2,7,^ = 1/2) + ^ I ^ = 1/2,7,^ = 1/2). 



Similarly, by omitting the label 7 and 7-^^- to simplify, 



tt"^ tt n) 



7r~ TT^ n) — 
I n" 7r° p) = 
I 7r°7r~p) = 



1 I 5/2,3/2) - I 3/2,3/2) - | 1/2,3/2) 



(12) 



'1 I 5/2,3/2) + ^^^13/2,3/2)- 



^ '1/2,3/2)-^ I 3/2,1/2) + ^ I 1/2,1/2) 



^15/2,3/2)+ ^ 



3V^ 



1 1/2,3/2)- 


1 
3 


3/2,1/2) + 


V2 
3 


1/2, 1/2) 


1 1/2,3/2)^ 


2 

3 


3/2,1/2)- 


x/2 
3 


1/2, 1/2) 
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Inverting the above equations we get, for example, 



I 1/2,1/2) = ^ I 7r° 7r%) - ^2 I 7r° TT" p) + ^2 I TT^ ttV) + 2 I TT" 7r+ n) j (13) 
|3/2,l/2) = ^(^v^ I 7r°7r°ra) - 2 I 7r°7r-p)- I 7r"7r» - I Tr^Tr+ra)^ 

In Fig. m we show the squared amplitude for J = 1/2 and /ttat =1/2 in the initial and 
the final state versus the total energy of the three-body system and the invariant mass of the 
meson-baryon subsystem formed by the second and third particle (ttN). A peak around an energy 
of 2100 MeV with a width of ~ 250 MeV appears when ^/s23 is close to 1670 MeV, thus having 
a 7riV*(1650) structure. The peak position and the width of this peak are compatible with the 
findings of various partial wave analyzes indicated by the PDG |28] about the A^*(2100), for which 
the peak position is found in the range 1855 - 2200 MeV and the width in the range of 69-360 
MeV. Thus we identify this peak with the A^*(2100). 

Since this peak appears when ^/s23 is close to the mass of the A^* (1650) and has been obtained 
using as input that vriV t-matrix which contains the information on the A^*(1650), we conclude 
that the inclusion of the iV*(1650) in the ttN subsystem is essential to generate a resonance at 
2100 MeV. 

In this former study, we do not find evidence for any resonance in the isospin 3/2 configuration, 
but the situation is different when we introduce coupled channels, as we discuss below. 



2.2 Inclusion of the ttKH, ttKA and tttjN channels 

Next, we solve the Faddeev equations with fourteen coupled channels: n'^ii'^n, 7r°7r~p, 7r°/C+S~, 
Ti^K^YP, TT^K^A, ttV, vr+TT-n, tt+K^J:', n-n+n, n-n\ ti-R+TP, tt~K^J:+, tt~K+A and n-rjp. 
Again, we label them as particles 1, 2 and 3 in the order in which they are written above. As 
there are no data for KT, KTj, KA KA, etc., we use the model of [23] to calculate the 
corresponding amplitudes. The vrA^ interaction below threshold is determined using the same 
model as for KT, and KA and above the threshold we use the experimental results. 

We continue to study those amplitudes where the isospin of the subsystem of particle 2 and 3, 
i.e., pion nucleon (and its coupled channels), is 1/2 in the initial as well as the final state. In Fig. 
|5]we show the ttt^N amplitude for total isospin J = 1/2 for such a case. 

As shown in Figj5]we obtain a peak at an energy of 2080 MeV with a width of 54 MeV for a a/s23 
near 1570 MeV, which we identify with the A^*(2100) listed in the PDG [2B]- Comparison of the 
Figs, m and E] shows that the inclusion of the Trie's, uKA and TTrjN channels makes the resonance 
more pronounced (by an order of magnitude in the squared T^-matrix) and much narrower. These 
changes in the results can be easily understood with respect to the previous ones obtained with 
only five coupled channels by noticing that now the wave function of the resonance contains extra 
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Figure 5: The A^*(2100) in the ttttN system including 14 coupled channels. 



12 



components which have smaller phase space in the decay of the resonance. At the same time, the 
TTTxN component becomes smaller due to the normalization of the wave function and, hence, the 
decay into -k-kN is also reduced. 

On exploring other isospin configurations we find a peak in the t^KK amplitudes with total 
isospin 7=3/2 and with isospin Ij^k = 1/2 in the initial as well as in the final state. In order to 
get this amplitude we have used the relation 

I / = 3/2, Ikk = 1/2) = ^(v^ I 7r°i^°A)+ | vr-iT+A)), (14) 

which has been obtained by writing the nKK states in isospin base analogously to Eqs. (fT2l[T3l) . In 
Fig. [6] we show the squared (J = 3/2, J/^a = 1/2 | T^{y/s, ^/s23) \ I = 3/2, Ika = 1/2) amplitude 
. A peak is found at a total energy of ~ 2126 MeV with ~ 42 MeV of width. In this case, the 
invariant mass ^/s23, at which the peak appears, is around 1590 MeV. This peak can be identified 
with the A(1910) listed in |28j, whose position, given by different partial wave analyzes, ranges up 
to 2070 MeV and the width varies from 190-500 MeV. 



|T*r|2(10-^ MeV"6) 




Figure 6: The A(1910) in the nKA system including 14 coupled channels. 
Thus the introduction of the ttKH, tcKA and TTfjN channels, together with the inclusion of 
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the iV*(1650) in the tiN t-matrix, is important to get this resonance. One should note that we 
get smaller widths than the experimental ones. The vrA^ decay channels are not considered in our 
approach and they should contribute to increase the widths. Note that this can be done even with 
a small ttN component, as implicitly assumed here, since there is more phase space for decay into 
the nN channel (see |T] for more discussion). 

We do not find any evidence of the A(1750), which could indicate a different structure for this 
state that the one studied in this work. 

3 Exploring the A^/o and A^ag systems by taking Nrnr, NKK 
and A^7r?7 as coupled channels 

Until now, we have investigated possible resonant states in the vrvrA^ system and its coupled channels 
which have been obtained by adding a pion to pseudoscalar-baryon systems which couple strongly 
in J'^ = l/2~ and isospin 1/2 configuration, i.e., ttN, KT,, KA and r]N. The invariant mass of this 
pseudoscalar-baryon subsystem has been varied around that of the A^*(1535) and A^*(1650), hence, 
treating the three-body system as a vrA^* system with 1500 < Mn* < 1760 MeV, although within 
the three-body Faddeev equations. There are other configurations of this three-body system, like 
A^ao(980) and iV/o(980), which we have not discussed so far. 

In order to study such a system, we must take NKK, Nim and Nt^t] as coupled channels, such 
that the vrvr and KK subsystem dynamically generate the /o(980) and the ttt] subsystem along 
with KK generates the ao(980) resonance. In this way, we can study the iV/o(980) and A^ao(980) 
systems simultaneously. Concretely, we take the following coupled channels into account: nvr^vr'^, 
p7r°7r~, mr^r], nTr"*"??", mi^Ti^, ptt^tt^, Pt^^Vj nK^ K~^ , nK^K^, pK^K~ . We label the particles as 
1, 2, 3 in the order in which they are written above. This means that the subsystem of particles 
2 and 3 consists of two pseudoscalar mesons whose invariant mass, y^S23, is varied around 980 
MeV. With these channels we solve the Eqs. ([3]) in the same formalism which we have explained 
in the previous sections. In this case we find that the NKK amplitude is bigger in magnitude 
as compared to those of the other coupled channels. We thus make isospin combinations of the 
NKK channels, similarly to Eqs. f|T2l [T3|) and obtain the amplitude for total isospin 7 = 1/2 and 
the isospin of the KK system, Ikr^ equal to or 1. 

In the case of total isospin of the NKK system equal to 1/2 with the isospin of the KK 
subsystem equal to one the amplitude 

(/ = 1/2, Ikk = 1\ mV~s, v^) I / = 1/2, Ikk = 1), (15) 

shows a peak around 2080 MeV , with a width of 51 MeV (which we do not show here), which we 
relate as the A^ao(980) partner of the peaks shown in Figs. [Hand [51 Thus the peak corresponding 
to the A^*(2100) has been seen in vrvriV system as well as in the NKK system. 
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Interestingly, along with this A^*(2100) state, we find another peak with even larger magnitude 
of the squared three-body amplitude at i/i = 1924 MeV with a width of 20 MeV. We show this 
peak in Fig. [7] for the NKK channel. 




Figure 7: A possible A^*(1910) in the NKK channels. 

This state is about 7 MeV below the NKK threshold (assuming an average mass for the kaons of 
496 MeV and 939 MeV for the nucleon). Therefore, this result indicates that the A^ao(980) system 
gets bound at around 1920 MeV. This possibility has been already suggested by the authors in [29] , 
in which they study the NKK channel using effective two-body potentials to describe the KN , 
KK, KN interactions. They find that the NKK system can get bound while the KK subsystem 
acts like the Oq. Our result is, thus, in agreement with the suggestions in [2^. Interestingly, the 
existence of a 1/2"*" A^* resonance around 1935 MeV has also been proposed earlier [30] on the 
basis of a study of the data on the 7p — >■ K~^A reaction in an isobar model, although other theories 
[3lj which include explicitly resonances up to 1855 MeV can reproduce these data (though further 
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work along these lines to include higher mass resonances is under way [32]). 

Since this peak found at 1920 MeV is below the three-body threshold and, in the two-body 
problem, the poles for the /o(980) and ao(980) appear below the KK threshold, the three particles 
in the system have associated complex momenta in the momentum representation. To avoid the use 
of unphysical complex momenta in the three-body system, which will lead to imaginary energies 
in the real plane, we give a minimum value, around 50 MeV, to the momentum of the particles. 
We have check the sensitivity of our results to the mentioned choice by changing the minimum 
momentum from 50 MeV to 100 MeV and we find the peak and width to remain almost unchanged. 

We also made the total isospin 1/2 combination of the NKK system by considering the KK 
subsystem in isospin in the initial and the final state, i.e., considering the A^/o(980) component 
of the NKK channel. In this case too, just as in the amplitude for A^ao(980), we find a peak 
around 1923 MeV with a width of 30 MeV and another one around 2052 MeV with a width of 60 
MeV. The magnitude of the A^/o(980) amplitude around 1920 MeV is very similar to the one of 
the iVao(980) amplitude (i.e., the one shown in Fig. [7]), but the magnitude of A^/o(980) amplitude 
around 2050 MeV is bigger than the magnitude of the A^ao(980) amplitude. 

From the whole study we would conclude that there are two iV*'s with J'^ = 1/2"*" in the energy 
region 1800 < < 2200 MeV. 

The peaks obtained in this approach are very neat and we associate them to physical resonance 
states. In the two body scattering it is customary to look for poles in the second Riemann sheet to 
associate them to resonances. The issue of poles for the three body problem within our formalism 
was addressed in the section VI of [5] . The difficulty to work with two complex variables, -y/s and 
\fs22, which induce complex three momenta needed in the evaluation of integrals are obvious. Yet, 
in [5] an approximate method was devised, which is suitable for the present context too, for the case 
when a subsystem of two particles can be treated as a resonance. Therefore the three-body system 
can be interpreted as a system of a particle and a resonance. This was the case for (/)/o(980) in [5] 
and the systems where resonances have been found in the present work can be treated similarly, 
for example, the NKK system where a resonance around 1920 MeV is found can be treated as a 
A^ao(980) system. In such cases, the problem can be reduced to a two-body scattering and usual 
poles can be identified in the complex energy (^i) plane. In as much as the resonances found 
here follow an approximate Breit-Wigner shape, as in [5J, the poles in the second Riemann sheet 
are guaranteed, as was shown in [5j. The amphtudes obtained in the present case indeed follow 
approximate Breit-Wigner distributions, and as in [5], they should have the corresponding poles 
in the y/s complex plane. 

4 Conclusions 

To summarize, we have extended our previous study of the vrvriV system and coupled channels [4], 
where the generation of the A^*(1710) was found, to higher energies. In this work the new input 
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is the experimental data on the vrA^ interaction where the information on excitation of both the 
A^*(1535) and the A^*(1650) is present, the latter of which was absent in our previous work [4J. Here, 
apart from confirming the A^*(1710), we find evidence for the other 1/2"'' A^*, i.e., the A^*(2100), 
and also for the 1/2"*" A(1910) resonance. The findings reported here indicate that the inclusion 
of the A^*(1650) in the interaction of the ttN subsystem is essential to generate these higher mass 
1/2"'" resonances. We have first made a search taking only the ttttN channels where a resonance 
having the properties of A^*(2100) was found. Later we included the vrKS, vrKA and vrr^S channels 
where the same resonance is produced but with larger magnitude and narrower width, indicating 
the addition of more channels to which the resonance couples strongly. No isospin 3/2 resonances 
is found in the study of the ttttN channels alone. However, the A(1910) is found on inclusion of the 
Tr/rS, ttKK and ttt^S channels. Further, we have investigated the NKK, Nirn and Nnt] channels 
where the KK — tttt subsystem rearranges itself as a /o(980) resonance, while KK — tti] acts like 
the ao(980). We obtain a new peak at ~ 1924 MeV, apart from the one corresponding to the 
A^*(2100), with a strong coupling to A^ao(980) and A^/o(980). Finally, we conclude this work by 
stating that the study of three-body systems, for the cases where a complete theoretical two-body 
input is not available, is also possible in our formalism using on shell experimental amplitudes. 
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